True ICC vs num
pooled ICCs
perf_crit |>
ggplot(aes(x = variance, y = rel_bias_icc_median, fill = model)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(var_combo_graph_f ~ icc_est_n,
labeller = label_bquote(rows = rho ==.(as.character(var_combo_graph_f)),
cols = k ==.(icc_est_n))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

ggsave("RPB_numICC.png", width = 6.8, height = 5, units = "in")
perf_crit |>
ggplot(aes(x = variance, y = rel_bias_icc_alt, fill = model)) +
geom_hline(yintercept = 1, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(var_combo_graph_f ~ icc_est_n,
labeller = label_bquote(rows = rho ==.(as.character(var_combo_graph_f)),
cols = k ==.(icc_est_n))) +
labs(x = NULL, y = "RPB - ratio", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)
ggsave("RPB_ratio_numICC.png", width = 6.8, height = 5, units = "in")
# perf_crit_MCSE <- perf_crit |>
# group_by(var_combo_graph, icc_est_n) |>
# summarise(mcse = max(bias_icc_mcse))
perf_crit |>
ggplot(aes(x = variance, y = bias_icc_median, fill = model)) +
geom_hline(yintercept = 0, linetype = "solid", colour = "black") +
geom_boxplot(alpha = .5) +
facet_grid(var_combo_graph_f ~ icc_est_n,
labeller = label_bquote(rows = rho ==.(as.character(var_combo_graph_f)),
cols = k ==.(icc_est_n))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

ggsave("PB_numICC.png", width = 6.8, height = 5, units = "in")
# ggplot(perf_crit, aes(var_combo_graph_f, bias_icc, col=model ) ) +
# facet_grid(variance ~ icc_est_n,
# labeller = label_bquote(rows = .(variance),
# cols = k ==.(icc_est_n))) +
# geom_point( alpha=0.10,
# position = position_dodge(width = 0.04) ) +
# geom_smooth( se=FALSE ) +
# geom_hline( yintercept = 0 ) +
# theme_minimal()
# ggplot(perf_crit, aes(var_combo_graph_f, bias_icc, col=skewness ) ) +
# facet_grid(variance ~ icc_est_n,
# labeller = label_bquote(rows = .(variance),
# cols = k ==.(icc_est_n))) +
# geom_point( alpha=0.10,
# position = position_dodge(width = 0.04) ) +
# geom_smooth( se=FALSE ) +
# geom_hline( yintercept = 0 ) +
# theme_minimal()
perf_crit |>
ggplot(aes(x = variance, y = RMSE, fill = model)) +
geom_hline(yintercept = .00, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(var_combo_graph_f ~ icc_est_n,
labeller = label_bquote(rows = rho ==.(as.character(var_combo_graph_f)),
cols = k ==.(icc_est_n))) +
labs(x = NULL, y = "RMSE", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

ggsave("RMSE_numICC.png", width = 6.8, height = 5, units = "in")
number clusters vs
cluster size
This is how the variance formulae perform across all generating ICC
values.
perf_crit |>
ggplot(aes(x = variance, y = rel_bias_icc, fill = model)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(nj_size ~ n_bar_size,
labeller = label_bquote(rows = j ==.(nj_size),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

#ggsave("RPB_j_nj.png", width = 6.8, height = 5, units = "in")
perf_crit |>
ggplot(aes(x = variance, y = bias_icc, fill = model)) +
geom_hline(yintercept = 0, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(nj_size ~ n_bar_size,
labeller = label_bquote(rows = j ==.(nj_size),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

#ggsave("PB_j_nj.png", width = 6.8, height = 5, units = "in")
perf_crit |> ggplot(aes(x = variance, y = RMSE, fill = model)) +
geom_hline(yintercept = .00, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(nj_size ~ n_bar_size,
labeller = label_bquote(rows = j ==.(nj_size),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RMSE", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

#ggsave("RMSE_j_nj.png", width = 6.8, height = 5, units = "in")
perf_crit |>
ggplot(aes(x = variance, y = rel_bias_icc, fill = nj_size)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(var_combo_graph_f ~ n_bar_size,
labeller = label_bquote(rows = rho ==.(as.character(var_combo_graph_f)),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

ggsave("RPB_j_nj_nomodel.png", width = 6.8, height = 5, units = "in")
perf_crit |>
ggplot(aes(x = variance, y = bias_icc, fill = nj_size)) +
geom_hline(yintercept = 0, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(var_combo_graph_f ~ n_bar_size,
labeller = label_bquote(rows = rho ==.(as.character(var_combo_graph_f)),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

ggsave("PB_j_nj_nomodel.png", width = 6.8, height = 5, units = "in")
perf_crit |> ggplot(aes(x = variance, y = RMSE, fill = nj_size)) +
geom_hline(yintercept = .00, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(var_combo_graph_f ~ n_bar_size,
labeller = label_bquote(rows = rho ==.(as.character(var_combo_graph_f)),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RMSE", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

#ggsave("RMSE_j_nj.png", width = 6.8, height = 5, units = "in")
perf_crit |>
ggplot(aes(x = variance, y = rel_bias_icc, fill = model)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(n_bar_prop ~ n_bar_size,
labeller = label_bquote(rows = zeta ==.(n_bar_prop),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(
legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |>
ggplot(aes(x = variance, y = bias_icc, fill = model)) +
geom_hline(yintercept = 0, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(n_bar_prop ~ n_bar_size,
labeller = label_bquote(rows = zeta ==.(n_bar_prop),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

filtered \(\rho = 0.50\)
perf_crit |> filter(var_combo_graph == 0.50) |>
ggplot(aes(x = variance, y = rel_bias_icc, fill = model)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(nj_size ~ n_bar_size,
labeller = label_bquote(rows = j ==.(nj_size),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(
legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |> filter(var_combo_graph == 0.50) |>
ggplot(aes(x = variance, y = bias_icc, fill = model)) +
geom_hline(yintercept = 0, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(nj_size ~ n_bar_size,
labeller = label_bquote(rows = j ==.(nj_size),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |> filter(var_combo_graph == 0.50) |>
ggplot(aes(x = variance, y = rel_bias_icc, fill = model)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(n_bar_prop ~ n_bar_size,
labeller = label_bquote(rows = zeta ==.(n_bar_prop),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |> filter(var_combo_graph == 0.50) |>
ggplot(aes(x = variance, y = bias_icc, fill = model)) +
geom_hline(yintercept = 0, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(n_bar_prop ~ n_bar_size,
labeller = label_bquote(rows = zeta ==.(n_bar_prop),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

filtered \(\rho = 0.90\)
perf_crit |> filter(var_combo_graph == 0.90) |>
ggplot(aes(x = variance, y = rel_bias_icc, fill = model)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(nj_size ~ n_bar_size,
labeller = label_bquote(rows = j ==.(nj_size),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(
legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |> filter(var_combo_graph == 0.90) |>
ggplot(aes(x = variance, y = bias_icc, fill = model)) +
geom_hline(yintercept = 0, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(nj_size ~ n_bar_size,
labeller = label_bquote(rows = j ==.(nj_size),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |> filter(var_combo_graph == 0.90) |>
ggplot(aes(x = variance, y = rel_bias_icc, fill = model)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(n_bar_prop ~ n_bar_size,
labeller = label_bquote(rows = zeta ==.(n_bar_prop),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |> filter(var_combo_graph == 0.90) |>
ggplot(aes(x = variance, y = bias_icc, fill = model)) +
geom_hline(yintercept = 0, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(n_bar_prop ~ n_bar_size,
labeller = label_bquote(rows = zeta ==.(n_bar_prop),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

filtered \(\rho = 0.05\)
perf_crit |> filter(var_combo_graph == 0.05) |>
ggplot(aes(x = variance, y = rel_bias_icc, fill = model)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(nj_size ~ n_bar_size,
labeller = label_bquote(rows = j ==.(nj_size),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(
legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |> filter(var_combo_graph == 0.05) |>
ggplot(aes(x = variance, y = bias_icc, fill = model)) +
geom_hline(yintercept = 0, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(nj_size ~ n_bar_size,
labeller = label_bquote(rows = j ==.(nj_size),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |> filter(var_combo_graph == 0.05) |>
ggplot(aes(x = variance, y = rel_bias_icc, fill = nj_size)) +
geom_hline(yintercept = .05, linetype = "dashed") +
geom_hline(yintercept = -.05, linetype = "dashed") +
geom_boxplot(alpha = .5) +
facet_grid(n_bar_prop ~ n_bar_size,
labeller = label_bquote(rows = zeta ==.(n_bar_prop),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "RPB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)

perf_crit |> filter(var_combo_graph == 0.05) |>
ggplot(aes(x = variance, y = bias_icc, fill = nj_size)) +
geom_hline(yintercept = 0, linetype = "solid") +
geom_boxplot(alpha = .5) +
facet_grid(n_bar_prop ~ n_bar_size,
labeller = label_bquote(rows = zeta ==.(n_bar_prop),
cols = bar(n_j)==.(n_bar_size))) +
labs(x = NULL, y = "PB", fill = NULL) +
theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(size = 6.5),
axis.text.y = element_text(size = 6.5),
plot.caption=element_text(hjust = 0, size = 6),
# legend.title=element_text(size = 6),
legend.text=element_text(size = 6.5),
axis.title.x = element_text(size=6.5),
axis.title.y = element_text(size=6.5),
# strip.background =element_rect(fill="white")
strip.text.x = element_text(size = 6.5, face="bold"),
strip.text.y = element_text(size = 6.5, face="bold")
)
